Local density of states in superconductor-strong ferromagnet structures 
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We study the dependence of the local density of states (LDOS) on coordinates for a 
superconductor- ferromagnet (S/F) bilayer and a S/F/S structure assuming that the exchange energy 
h in the ferromagnet is sufficiently large: hr » 1, where r is the elastic relaxation time. This limit 
cannot be described by the Usadel equation and we solve the more general Eilenberger equation. We 
demonstrate that, in the main approximation in the parameter (hr) -1 , the proximity effect does not 
lead to a modification of the LDOS in the S/F system and a non-trivial dependence on coordinates 
shows up in next orders in (/it) -1 . In the S/F/S sandwich the correction to the LDOS is nonzero 
in the main approximation and depends on the phase difference between the superconductors. We 
also calculate the superconducting critical temperature T c for the bilayered system and show that 
it does not depend on the exchange energy of the ferromagnet in the limit of large h and a thick F 



?H ' layer. 
03 \ I. INTRODUCTION 



In recent years the interest in the proximity effect between superconducting (S) and ferromagnetic (F) layers has 
increased considerably. The actual progress in the preparation of high quality metallic multilayer systems allows a 
careful study of the mutual interaction of superconductivity andrferromagnetism in hybrid S/F structures. The prox- 
q ' imity effect manifests itself in e.g. changes of the density of statesu (DQSJ or in the dependence of the superconducting 
critical temperature T c on the thickness of the F-layer (see e.g. Refsflcl). 

Many theoretical works have been devoted to the study of such structures, which is not a simple task. In most of 
them the DOS and the critical temperature of the superconducting transition T c were analyzed in the "dirty limit" . 
■^j- ■ This means that the retarded and advanced Gteea's functions, which determine the thermodynamical properties, were 
P\| ' obtained from the Usadel equation (e.g RefsHTlJ) which is simpler than the more general Eilenberger equation. In 
other words, it was assumed that the mean free path I is much shorter than any characteristic length (with exception 
of the Fermi wave length) , and that all energies involved in the problem are smaller than r 1 , where r is the elastic 
relaxation time. In particular, the condition hr -C 1 has to be satisfied (h is the exchange energy). It was shown that 
in this limit the DOS in the F layer oscillates and decays with increasing exchange energy h. In the limit hr <C 1, the 
period of the oscillations and the decay length are comparable to each other apd-are of the order y/D/h. ■— ■ 

Calculations in the opposite (purely ballistic) case were performed in RefsOtJ. Haltermanns and VallsEa solved 
the Bogolyubov and selfpCnnsistency equations numerically and presented the results for the LDOS in a F/S structure. 
I i Zareyan and co-workersEil calculated the LDOS of a F film in contact with a superconductor. They assumed that 
the electrons in the F film are scattered only at the rough F/Vac surface (Vac stands for Vacuum). However, not so 
much attention has been payed to the more realistic case hr 1. Except for Ref.Ej, where this limit was considered 
for very thin ferromagnetic films, a study of the proximity effect in disordered ferromagnetic conductors with strong 
exchange energies h is still lacking. This limit is very important when dealing with ferromagnetic layers made of 
transition metals, as Fe or Ni, for which h < leV, and whose sizes are larger than the mean free path /. 

In the present paper, we analyze the DOS and the critical temperature T c of the superconducting transition in S/F 
structures for which the condition 



hr > 1 (1) 

is satisfied. Here r is the elastic relaxation time due to impurity scattering in the bulk. According to condition (|l|) 
the exchange interaction, which is proportional to h is much stronger than the interaction with impurities, which is of 
the order of r _1 . We will show that in this limit the decay length of the superconducting condensate function induced 
in the ferromagnet is of the order of the mean free path I. Therefore, despite its smallness, it is very important to 
retain the impurity scattering term if the thickness of the F layer is larger than On the other hand, the product 
tT c may be smaller or larger than unity. We will pee-that in the limit (|l|) thermedmpamical quantities, such as DOS 
and T c , differ considerably from those in the dirtyB~li3 and purely ballistic limitaiitHl. 
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As it has been pointed out in Refs.t-fcl, the Usadel equation is not applicable in the limit ([!]), and the more 
general Eilenberger equation should be solved. One can solve the latter equation with the help of an expansion in the 
parameter (hr)^ 1 for both the strong and weak proximity effect. In the case of a weak proximity effect the solution 
can be obtained even for an arbitrary impurity concentrationEj, i.e. for an arbitrary value of the parameter hr. It will 
be shown that in the limit of large hr the condensate function / oscillates in space with a period vp/h and decays 
on a [length of the order of the mean free path /. If one neglects the impurity scattering in the bulk, as was done in 
RefslIJO, the results for DOS depend essentially on the boundary conditions at the F/Vac interface. This is because 
in that case the condensate function / is formed by interfering waves reflected from the F/Yac boundary. .Therefore, 
the thermodynamical properties depend sensitively on whether these reflections are specularlia or diffusivetll. 

In the present paper, we assume that the thickness of the ferromagnetic layer dp is much greater than the mean 
free path I and therefore reflected waves are unimportant provided the bulk impurity scattering dominates. We show 
that in the main approximation in the parameter (/it) -1 , the LDOS in the ferromagnetic region in a F/S bilayer is 
not affected by the proximity effect. On the other hand the amplitude of the condensate function in the ferromagnet 
near the S/F interface is equal to that in the superconductor S if the interface transmittancc is high, and decays away 
from the interface at distances of the order I. In the limit (|l|), the mean free path I is much larger than the period of 
the oscillations (~ Vp/h) and may be comparable with dp. Therefore, one can speak about long-range penetration 
of the superconducting condensate into the ferromagnet. In the presence of two superconductors (e.g. in a S/F/S 
sandwich) the situation changes: a correction to the DOS in F arises in the main approximation and its magnitude 
depends on the phase difference between the superconductors and on the thickness dp of the F layer vanishing in the 
limit dp — > oo. We also analyze a variation of the critical temperature of the superconducting transition T c for a 
S/F bilayer for which the condition ([j]) is satisfied. We find the dependence of T c on the thickness d$ of the 
superconductor in the limit dp 3> I and compare it with the results obtained in the dirty limit. 



II. BASIC EQUATIONS 

Physical quantities of interest can be computed using quasiclassical retarded g R and advanced g A Green's functions. 
In order to find the retarded Green's function g R (for brevity we omit the index R) one should solve the Eilenberger 
equation (the advanced Green's function can be found in the same way) 

fiv F d x g - i(e + h) [f 3 ,g] - i \L,g] + (1/2t) [(g), g] = . (2) 
Here g is the retarded Green function in the Nambu space. It obeys the normalization condition 

9 2 = 1 • (3) 

The angle brackets denote the averaging over angles: ((...)) = f Q d(j,(...), /i = cos# , 9 is the angle between the 
momentum and the x-axis (in all the S/F systems considered below the ir-axis is perpendicular to the S/F interfaces), 
and vp is the Fermi velocity. The exchange field h vanishes in the S-layers. A is the matrix 



A = 



A 

-A* 



and A is the pair potential in the superconductor. We assume that in the F-layers the electron-electron interaction 
vanishes and therefore A = 0. The last term of Eq.(|^) describes the. effect of nonmagnetic impurities. Eq.(j|) is 
supplemented by the Zaitsev boundary condition at the S/F interfaced: 

&\R- Ra 2 + - .s 2 ) 2 ) = j[h, h] ■ (4) 

Here a and s are the antisymmetric and the symmetric (in /j) parts of the Green's function g; R and T are the 
reflection and transmition coefficients. Eqs. (0fl) describe the system completely and are the basic equations in our 
study. In spite of a rather simple symbolic representation they are quite complicated and a general solution can hardly 
be obtained for an arbitrary impurity concentration. Accordingly, further calculations are performed under certain 
assumptions that allow us to simplify equations (0-0) for values of parameters which still correspond to real systems. 
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III. LOCAL DENSITY OF STATES IN A S/F BILAYER 



In this section, it is assumed that the mean free path l s in the superconductor is larger than the coherence length 
£ s . In this case the last term in Eq. (||) can be disregarded in the S region. We consider first a S/F bilayer. The S/F 
interface is located at x = 0, the superconductor and ferromagnet occupy the regions x < and x > 0, respectively. 
We also assume that the thickness of both layers are much larger than the characteristic lengths over which g varies, 
i.e £ s in the S-layer and I in the F layer. We focus first on the LDOS in the superconductor for subgap energies: e < A. 
It is clear that in this energy range and for \x\ ^> £ s the LDOS vanishes as it should be in a bulk superconductor. 

We assume a perfect S/F interface transparency. In this case and according to the boundary condition Eq. (^) 
both the symmetric and the antisymmetric parts must be continuous at the interface. For simplicity we approximate 
the dependence A(a;) by a step like function A(x) = AQ(—x). Then, the solution of Eq. (||) in the superconducting 
region has the form 

9S = (e/£ - sgn M • A/£Ce K «*) f 3 + (A/£ - sgn/i • (e/0Ce^ x ) if 2 + Csgnfie^h . (5) 



Here £ = \J t 2 ~ A 2 , n s — —2i^/\fi\vp and C is a constant which has to be determined from the boundary conditions. 
Note that the solution (|J) can also be represented in terms of the eigenfunctions U = {u, v) of the Bogolyubov 
equation, i.e. it can be written as 

9S = gsb + f 3 U + {x)® U(x) , (6) 

where g s t> is the Green function in the bulk. For energies |e| < A the functions U(x) decay exponentially away from 
the S/F interface. 

Now we have to solve Eq. (U) in the ferromagnet. In the limit determined by Eq. ([[]) one can check that the 
solution of Eq. (0) in the main approximation with respect to (hr)^ 1 is 

Qf = T3 + .Dexp(— nx/\fji\l)iT2 + sgn/i • Dexp(— Kx/\fj,\l)fi , (7) 

where n = 1 — 2i(e + h)r. The second and third terms in Eq. (0) describe the condensate function induced in the 
ferromagnet due to the proximity effect. 

As opposed tcj-piany authors claiming that the decay length of the condensate function is proportional to Vp/h (see 
for example Ref.Q), Eq. (Q) clearly shows that the penetration depth of the condensate is of the order of the mean free 
path I and does not depend on the exchange field h. Thus, for clean and strong ferromagnets, the superconducting 
condensate penetrates over long distances compared to the magnetic length vp/h. As mentioned in the introduction, 
in the dirty case, hr <C 1, the penetration length of the condensate into the ferromagnet decreases with increasing 
h, which is true for the singlet component of the condensate function Ilha triplet component of the condensate may 
penetrate into the ferromagnet over a long length of the order yj Z?/TtL3'ti This triplet component may be induced 
by the presence of a magnetic inhomogeneity close to the S/F interfacai3) . However, we see from Eq. (0) that a 
long-range proximity effect exists even for singlet pairing and homogeneous magnetization. Using the fact that g is 
continuous at x = 0, one can determine the constants D and C from Eqs. (||) and (0) 

C = D= Zl±f (8) 
A v ; 

The normalized LDOS is given by the well known expression 

v = v/u Q = (l/2)ReTrf 3 g , (9) 
where z/ is the LDOS in the normal state. From Eqs. &0) one obtains in the S region (e < A) 



v{x,e) =<exp(2 v / A 2 - e 2 x/\n\v F ) > . (10) 

In Fig.^, we plot the spatial dependence of v and the absolute value of the imaginary part of the symmetric condensate 
function (second term in Eq. (0)). Due to the proximity effect of the ferromagnet the superconductor becomes gapless 
in the region close to the S/F interface. According to Eqs. (0) and (||) the LDOS in the ferromagnet remains unchanged 
despite the presence of the condensate function induced in the F region (second and third term in Eq. (0)). The 
reason for that is the cancellation of the symmetric (second term in Eq. (ffl)) and antisymmetric parts (third term in 
Eq. (Q)) in the normalization condition (||). Thus, the coefficient in front of T3 in Eq. (0) is 1 and therefore v = 1. 
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FIG. 1. The spatial dependence of the normalized LDOS v (dashed line), and the imaginary part of the symmetric condensate 
function |Im/2| (solid line). The S (F) layer is located in the region x < (x > 0). Ao is the value of A for T = 0. A/Ao = 0.8, 
e/Ao = 0.4,rAo = 1 and h/ Ao = 10. Although the function |Im/2| exhibits a strongly oscillating and weakly decaying behavior 
in the ferromagnet, the LDOS remains unchanged in this region. 



The correction to the LDOS in the ferromagnet is not zero if one takes into account higher order terms in the 
expansion in the parameter (/it) -1 . We calculate this correction assuming for simplicity that the transparency of the 
S/F interface is low. In this case, the Green function in the superconductor is not affected by the proximity effect 
and is given by 

9s = G s f 3 + F s iT 2 , (11) 

where Gs = e/Ve 2 — A 2 and F$ = A/\/e 2 — A 2 . For a low interface transparency, the boundary condition, Eq. (Q), 
reduces to 

a = 7 [sf, ss] — iFsh . (12) 

Here 7 = T (fi) / 4R(fi) <Cri.is the parameter describing the transmittance of the interface. The way how to solve Eq. 
(||) was presented in Ref.li3. The exact Green function is given by 

9f(x) = f 3 + h{x)%Vi + fi(x)n = f 3 + J (dfc/27r)(/ 2fe if 2 - (jjl / 1 K)hkTxd x )e lkx , (13) 

where 

hk = M[1 ^L/A/)] hK (w(1/M) ~ (^/ M » + 7M] , (14) 

M = (klii) 2 + k 2 and k = 1 — 2i(e + h)r. In the case under consideration, i.e hr ^> 1, one can easily show that 
\k(1/M)\ < 1 and hence we can expand Eq. dlj): 



2nFsl 

h* = (7M + «7oA^=i) • (15) 



For simplicity we have assumed that the transmission coefficient T(/i) has a sharp maximum at 9 = 0, i.e. 
7(/x) = 7o<5(m — 1). The first correction to the normalized LDOS proportional to (/it) -1 can be now obtained. 
After cumbersome but straightforward calculations one obtains at low energies 
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= 2 Re ((/2) 2 - (A) 2 ) - -7o(e" 2a5/ 74M sin^rx/O 



(16) 



Eq. ( |l6|) shows that for sufficiently large hr, the correction to the LDOS of the normal metal is small. This cpsrection 
to the LDOS proportional to (/it) -1 reveals a damped-oscillatory behavior similar as the one reported in Ref.t2l, where 
the case hr <C 1 (dirty limit) was considered. The spatial dependence of the local DOS is shown in Fig.||. 




FIG. 2. Spatial dependence of the LDOS in the ferromagnet for different values of hr. The solid, dashed and dash-dotted 
lines correspond to hr = 15, hr = 10 and hr = 5, respectively. 



IV. LOCAL DENSITY OF STATES IN A S/F/S SANDWICH 

Now we consider a Josephson-like S/F/S structure (see Fig.|J). In this case, the correction to the DOS is not zero 
even in the main approximation in the parameter (hr)^ 1 . The thickness of the F-layer is dp , and we assume again 
that the transparency of the S /F interfaces is low. The solution of Eq. (Q) can be sought in the following form 

g = r 3 + s + a, (17) 

where s and a are the symmetric and antisymmetric part of the condensate function induced in the ferromagnet. The 
equations which determine s and a can be obtained from Eq. (^) and written in the form 

(nl) 2 dxJ- k 2 s- k(s) =0 
a = -(pl/n)t 3 d x s , 

where k is defined in Eq. (|l4|). 
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x=-d F /2 



x=d F /2 

FIG. 3. The S/F/S structure. 



In the limit of a weak proximity effect, the boundary conditions at x = ±d F /2 are [the upper (lower) sign corresponds 
to x = d F /2 (x = -dp/2)] 



s ■ 



(19) 



Fs is the condensate function of the superconductors which now is given by Fs(y>) — ir^Fs exp(±iipt 3 ) , where F$ is 
defined in Eq.([TT|) and ip is the phase difference between, the superconductors. The general solution of Eq. (18) with 
the boundary conditions Eq.(jl^) was presented in RefEj. In the limit hr ^> 1 one obtains 



i(x) = -2M/i7Fs^(i-r 2 cos^/2 + if 1 (-l) n sm^/2) / ^ e -ikx e ik(2n+i)d F j M 

n 

t\ ■ ~ ( ir\\ cosh 9(x) , ■ - • / /a\ sinh 0(x) 

= -jF s it 2 cos(^/2) s . nh Kd y /2) + tn srn(y/2) cosh e{ ^> 2) 



(20) 



where 6(x) — nx/l\fi\ and n takes integer values from -co to +oo. From Eq. (|18[) one obtains the expression for a 



a(x) = 7-Fs 



f 2 sin((^/2) 



cosh0(x) 
cosh(9(d F /2) 



fi cos(</?/2) 



sinh 0(a;) 



smh6(d F /2) 



The correction to the LDOS due to the proximity effect is then given by 



8v 



1 



Re(s 2 



A 2 



r Re( 
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sinh 2 29(d F /2) 



(1 + cos <p cosh 20(d F /2))) , 



(21) 



(22) 



where 28(d F /2) = Kd F As one could expect, 8v vanishes in the limit of a very large thickness d F of the F-layer. 
At the same time, one can see that the correction to the LDOS is not x-dependent but it depends on the phase 
difference ip. This behavior is quite interesting and rather unexpected. The correction to the LDOS, Eq. (p2|), may 
be both positive and negative depending on the phase (p. In Fig.|] and Fig ^5] we plot the dependence of 8v on the 
thickness d F and on the phase difference <p. As before we assume that 7 2 (/i) = (7o/<5/i)<5(/i— 1), where <5/i is the 
width of the peak in the dependence 7(/i). Thus, for a S/F/S sandwich the correction to the LDOS 5v 
is much larger than in the case of a S/F bilayer. In the latter case the correction is proportional to the 
small factor (hr)^ 1 , while in the S/F/S structure 5v is finite in zeroth order in (/it) -1 and proportional 
to the large factor 1/S/j,. 
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FIG. 4. LDOS versus the thickness dp of the F-layer. Here er = 0.1, At — 1 and tp — 0. The solid and dashed lines 
correspond to hr — 15 and hr = 10. 
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FIG. 5. LDOS versus the phase difference <p for the parameters er = 0.1, At — 1 and (If/1 = 5. The solid, dashed and 
dash-dotted lines correspond to hr — 15, hr = 10 and hr — 5, respectively. 



V. CRITICAL TEMPERATURE FOR A S/F BILAYER 



Another thermodynamic quantity we are interested in, is the supercymducting critical temperature of a S/F bilayer 
in the limit hr 1. The opposite case (hr <C 1) was studied in RefsJj Here we consider a bilayer system consisting 
of one S-layer of thickness ds and one F-layer of thickness dp. Of course, the superconducting transition temperature 
can considerably change only provided the thickness of the superconducting layer is not very large d s < £s . We assume 
that the S/F interface has a very high transparency, otherwise the proximity effect does not affect substantially the 
critical temperature T c . The computation of the condensate function / = s + a from Eq. (^) follows as before. In 
the limit hr 1 the solution in the F-layer is given simply by Eq. ([?]). The solution in the S region requires a little 
more care. For temperatures close to T c one can seek for solutions in the form 

9S = f 3 + F 2 if 2 +Fin . (23) 

If we assume that A (a;) is a slow varying function of x, then the function F 2 can be written as0 

F 2 = A/e + SF (24) 
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Bearing-in mind that the functions F\ and F2 are continuous at the S/F interface, one obtains the following equation 
for SFB 



(niydiJF - n s 5F - k s (SF) = -2\{i\Ik s D8{x) , (25) 

where k s = — 1 + 2ier, and D is the constant from Eq. (^). Equation ( p5| ) can be easily solved in the Fourier space. 
The solution is given by 

91k D 

M s (1 + K s (l/M s )) 

Here Ms = (kl/i) 2 + n 2 s . In the dirty limit, i.e. when er -C 1, the expression Eq. ( p6| ) can be simplified 

3ln s D 

SFk = W ~ 6»er ■ ( } 

and the value of the function SF(x) at x = is given by 



From Eqs. (Q), and the fact that all the functions are continuous at x — 0, one can determine the constant 

D 

D=* £E . (29) 

We see that the constant D does not depend on the strength of the exchange field h, and therefore the condensate 
function in the superconductor is not /i-dependent. Thus, we come to the following conclusion: in the limit /it > 1 
and if (if > I the thermodynamic quantities of the superconductor do not depend on the strength of the exchange 
field h; in particular, the critical temperature of the bilayer does not depend on the thickness of the ferromagnetic 
layer. This result is in qualitatively agreement with experimental data presented in Ref.p, in which the dependence 
T*(d,Gd) of Nb/Gd systems was determined and a saturation of T* for large dad was observed. 

In order to estimate the critical temperature T* of the system, we need to calculate the pair potential A (a;). The 
latter satisfies the self-consistency equation 

A = \ J detanh(e/3)(F.f - F*) (30) 

We notice that all equations presented above are for the retarded component. The advanced Green function F£ can 
be obtained using the relation (F^'Y — ~ ^2 ■ O ne can easily check that at x = the integral in Eq. ( |30| ) diverges as 
1/y/t and therefore the function A(x) must vanish at x — 0. Thus, the boundary condition for the pair potential at 
the S/F interface can be written as 

A(0) = . (31) 

Using the Ginsburg-Landau equation 

e GL A" = -A , (32) 



where ^gl = £o/v 1 — T*/T c , and the boundary condition Eq. ( |3l| ) we obtain 

A = Asm{x + d s )/£ GL . (33) 
At the boundary with the vacuum (x — —d s ) d x A — 0. From this condition we obtain an expression which determines 

1 C 

1 - T*JT C = « /24) 2 • (34) 

We remind that equation ( |34| ) is only valid in the case \T* — T c \ < T c . It gives the asymptotic value of T* in 
the limit hr ^> 1 and dp ^> I- This formula describes the usual proximity effect in a S/N structure: T* coincides 
with T c ii d s >> £0 an d decreases with increasing d s Ref.Eil. The same expression for T* is valid in the dirty 
limit (hr < 1, see Ref.Q). 
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VI. CONCLUSION 



In conclusion, we have considered different S/F structures in the case of a strong exchange field h, whenJthii 
condition hr ^> f is satisfied,— The results of our analysis differ considerably from those obtained in the "dirty"ETE£l 
and in the pure ballistic limited. We have shown that the thermodynamic properties ( such as local DOS and critical 
temperature) of a S/F bilayer do not depend on h in the main approximation with respect to the small parameter 
(/it) -1 . It is worth mentioning that in this approximation the LDOS does not change in the ferromagnet below T c . 
The superconducting gap is however suppressed near the S/F interface (gapless state) and restored at distances of 
the order of £ s ~ vp/A. In a S/F/S sandwich the correction to the LDOS in the ferromagnet is nonzero, spatially 
constant, and depends on the phase difference tp between the superconductors. The different behavior of the LDOS 
in a S/F and S/F/S structure is due to the interference of the induced condensate functions created at each S/F 
interface. Thus, the changes of the LDOS v due to the proximity effect may be observed more clearly in a S/F/S 
structure by measuring the dependence of v on the phase difference (p between the superconductors. 

We would like to thank SFB 491 Magnetische Heterostrukturen for financial support. 
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